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Abstract 

The classical mechanics, exact quantum mechanics and semiclassical quantum mechanics of the 
billiard in the triaxial ellipsoid is investigated. The system is separable in ellipsoidal coordinates. 
A smooth description of the motion is given in terms of a geodesic flow on a solid torus, which 
is a fourfold cover of the interior of the ellipsoid. Two crossing separatrices lead to four generic 
types of motion. The action variables of the system are integrals of a single Abelian differential of 
second kind on a hyperelliptic curve of genus 2. The classical separability carries over to quantum 
mechanics giving two versions of generalized Lame equations according to the two sets of classical 
coordinates. The quantum eigenvalues define a lattice when transformed to classical action space. 
Away from the separatrix surfaces the lattice is given by EBK quantization rules for the four types 
of classical motion. The transition between the four lattices is described by a uniform semiclassical 
quantization scheme based on a WKB ansatz. The tunneling between tori is given by penetration 
integrals which again are integrals of the same Abelian differential that gives the classical action 
variables. It turns out that the quantum mechanics of ellipsoidal billiards is semiclassically most 
elegantly explained by the investigation of its hyperelliptic curve and the real and purely imaginary 
periods of a single Abelian differential. 

PACS: 03.20.+i, 03.65.Ge, 03.65.Sq 
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1 Introduction 

Almost exactly 160 years ago, Carl Gustav Jacobi was able to separate the geodesic flow on 
ellipsoidal surfaces [Q]. In his letter from December 28, 1838 to his colleague Friedrich Wilhelm 
Bessel he wrote: 



li Ich habe vorgestern die geodatische Linie fur ein Ellipsoid mit drei ungleichen Achsen auf 
Quadraturen zuriickgefiihrt. Es sind die einfachsten Formeln von der Welt, Abelsche Integrale, die 
sich in die bekannten elliptischen verwandeln, wenn man zwei Achsen gleich setzt." 1 

The billiard motion inside an n-dimensional ellipsoid appears as the singular limit of the 
geodesic flow on an (n + l)-dimensional ellipsoidal surface with one semiaxis approaching zero. 
The starting point in Jacobi's treatment is what nowadays is called Hamilton-Jacobi ansatz. By 
introducing ellipsoidal coordinates Jacobi has shown that the integration of the Hamilton-Jacobi 
generating function leads to Abelian integrals. From Jacobi's point of view this insight essentially 
solves the separation problem. Jacobi's integrals represent the Abel transform of the ellipsoidal 
coordinates for which the time evolution is trivial. To give explicit expression for the time 
evolution of the ellipsoidal coordinates themselves it is necessary to invert the Abel map. The 
solution of this problem, the so-called Jacobi inversion problem, requires deep insight in the theory 
of meromorphic functions on hyperelliptic curves and has given rise to the definition of theta 
functions Q. This area constituted a highlight in 19th century mathematics. With the advent of 
quantum mechanics the attention of the scientific community was shifted from these non-linear 
finite dimensional problems to linear infinite dimensional problems. Recently classical mechanics 
has experienced a revival with two main directions. On the one hand computers have induced 
a boom in the study of non-integrable systems, essentially by allowing for the visualization of 
chaotic phenomena like the break up of Kolmogorov-Arnold-Moser tori. The quantum mechanics 
of non-integrable systems today is a main topic in physics. On the other hand the investigation 
of soliton equations has given deep insights into the theory of integrable systems and a lot of the 
knowledge about integrable systems of the 19th century has been revived. 

In this paper on ellipsoidal quantum billiards we explain the quantum mechanics of an integrable 
system in terms of the corresponding classical system via a semiclassical approach. The main object 
will be the hyperelliptic curve of Jacobi's classical theory. As usual the curve comes into play in 
order to give a defintion of the action differential. The action differential corresponding to the 
ellipsoidal billiard defines a hyperelliptic curve of genus 2 on which it is an Abelian differential 
of second kind. The real and purely imaginary periods of this differential enter the semiclassical 
quantization scheme in a very natural way. The presentation of this unified view of classical and 
semiclassical treatment is the main theme of this paper. 

According to the Liouville-Arnold theorem the phase space of an integrable system with / 
degrees of freedom is foliated by invariant manifolds which (almost everywhere) have the topology 
of /-tori. The most elegant phase space coordinates are action-angle variables (I,<p), where the 
action variables i" label the tori and the angles <p parametrize the torus for fixed i". With the original 
phase space variables (p, q) the action variables I are obtained from integrating the Liouville 1-form 
p dq along / independent cycles 7, on the torus according to 

h = — I pdq, i = 1,...,/. (1) 



0, (2) 



Hamilton's equation reduce to 

. = dH(I) 

1 d(fi 

1 English translation: The day before yesterday, I reduced the geodesic line of an ellipsoid with three unequal axes 
to quadratures. The formulas are the simplest in the world, Abelian integrals, transforming into the known elliptical 
ones if two axes are made equal. 
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dH(I) . 

fi = qj = Ui , i = l,... J (3) 

with u>i the constant frequencies. The time evolution becomes trivial. Although the importance of 
action-angle variables is stressed in any text book on classical mechanics, especially as the starting 
point for the study of non-integrable perturbations of integrable systems there can be found only 
few non-trivial examples in the literature for which the action variables are explicitely calculated. 
P. H. Richter Q started to fill this gap for integrable tops and recently this presentation has been 
given for various systems, e.g. for the Kovalevskaya top H integrable billiards with and without 
potential |?], ||, H the integrable motion of a particle with respect to the Kerr metric p~Tf] , and 
the motion of a particle in the presence of two Newton potentials - the so-called two-center-problem. 
It turns out that the presentation of energy surfaces H(I) = E in action space may be considered 
as the most compact description of an integrable system Jl2[ . 

The importance of action variables extends to quantum mechanics in the following way. In a 
semiclassical sense the Liouville-Arnold tori carry the quantum mechanical wave functions. Sta- 
tionary solutions of Schrodinger's equation result from single valuedness conditions imposed on the 
wave functions carried on the tori. These give the semiclassical quantization conditions 

Ii = {m + ai/4:)h, i = l, ...,/, (4) 

with quantum numbers rij and Maslov indices Qj. The are purely classical indices of the cor- 
responding Liouville-Arnold torus which is a Lagrangian manifold. They come into play because 
quantum mechanics is considered with respect to only half the phase space variables (p, q) - usually 
in configuration space representation, i.e. with respect to the qi [|lj|. The Maslov indices character- 
ize the singularities of the projection of the Lagrangian manifold to configuration space which lead 



to phase shifts of semiclassical wave functions supported on the tori [14, |13[| . The Maslov indices 
depend on the choice of the cycles ji in Eq. <^j. In the case of a separable system and a canonical 
choice of the cycles on the torus according to 

ji:dqj = 0, j j^i, (5) 

we simply have oti = if the ith degree of freedom is of rotational type and Oj = 2 if the ith degree 
of freedom is of oscillatory type. The EBK quantization ([|) was the center of the old quantum 
mechanics of Bohr and Sommerfeld before 1926. The fact that this quantization assumes that 
phase space is foliated by invariant tori was realized by Einstein, but his 1917 paper on this matter 
fli~5l| was hardly recognized at that time. 

The phase space of the ellipsoidal billiard is foliated by four types of tori which have different 
Maslov indices. Two crossing separatrix surfaces divide the action space into four regions - one 
four each type of tori. This means that the simple EBK quantization condition (Q) is not uniformly 
applicable to the ellipsoidal billiard: the quantum mechanical tunneling between the different types 
of tori has to be taken into account. Both effects can semiclassically be incorporated by a WKB 
ansatz for the wave function. The tunneling between tori is then described by tunnel matrices which 
connect the amplitudes of WKB wave functions in different classically allowed configuration space 
areas. The main ingredient for the tunnel matrices is a penetration integral. For the ellipsoidal 
billiard there exist two such penetration integrals - one for each separatrix. 

The differentials for both penetration integrals are identical. They are even identical to the 
differential for the action integrals of the three degrees of freedom, the only difference is the in- 
tegration path. The action and penetration integrals therefore appear as the real and purely 
imaginary periods of a single Abelian differential of second kind. This is how semiclassical quan- 
tum mechanics extends the meaning of the originally classical hyperelliptic curve and how quantum 
mechanics appears as a "complexification" of classical mechanics. 

Within the last few years the study of billiards has become very popular in connection with 
the investigation of the quantum mechanics of classically chaotic systems. The quantum mechanics 
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of two-dimensional billiards can easily be investigated experimentally by flat microwave cavities 
for which one component of the electric field vector mimics the scalar quantum mechanical wave 



function [16, 17, 18]. The relation between Schrodinger's equation for a quantum billiard and 



Maxwell's equations for the electromagnetic field in a three-dimensional cavity is complicated by 



the vector character of the electromagnetic field [19]. Three-dimensional billiards have a direct 
physical interpretation as models for atomic nuclei [20] and metal clusters [21]. Recently their 
importance has been rediscovered in connection with lasing droplets [22|. The semiclassical analysis 
of rotationally symmetric ellipsoids can be found e.g. in [p3| , |24| , p5|j . 

This paper is organized as follows. In Section § we summarize the classical aspects of the 
ellipsoidal billiard. We introduce constants of the motion, discuss the different types of tori and give 
a regularization of the ellipsoidal coordinates. In Section || the hyperelliptic curve associated with 
the ellipsoidal billiard is investigated. The separated Schrodinger equation is solved in Section ^. 
In Section |5| a uniform semiclassical quantization scheme in terms of a WKB ansatz is performed 
and a representation of the quantum eigenvalues in classical action space is given. In Section |6| we 
comment on how the degenerate versions of the ellipsoidal billiard, i.e. the prolate and the oblate 
ellipsoidal billiard and the billiard in the sphere, appear as special cases of the general triaxial 
ellipsoidal billiard. We conclude with some brief remarks and an outlook in Section |7[ 

2 The Classical System 

We consider the free motion of a particle of unit mass inside the general triaxial ellipsoid in K 3 
defined by 

* 2 + Ti + li = i W 

with < b < a < 1. The particle is elastically reflected when it hits the boundary ellipsoid. 
Throughout this paper we take (a, b) = (0.7, 0.3) in our numerical calculations. 

Hamilton's equations of motion and the reflection condition are separable in ellipsoidal coordi- 
nates (£, rj, C)- Each of them parametrizes a family of confocal quadrics 

^ + -^- + -^- = 1, (7) 

s 2 s 2 — b 2 s 2 — a 2 

where s G For 1 > s = £ > a all terms in Eq. (0) are positive and the equation 

defines a family of confocal ellipsoids. Their intersections with the (x, y)-plane, the (x, z)-plane 
and the (y, z)-plane are planar ellipses with foci at (x,y) = (±6,0), (x,z) = (±a, 0) and (y,z) = 
(±(a 2 — 6 2 ) 1 / 2 ,0), respectively. For a > s = rj > b the third term in Eq. (|7]) becomes negative. 
Eq. @ thus gives confocal one sheeted hyperboloids. Their intersections with the (x,y)-plane are 
planar ellipses with foci (x, y) = (±6, 0); the intersections with the (x, z)-plane and the (y, z)-plane 
are planar hyperbolas with foci at (x,z) = (±o, 0) and (y,z) = (±(a 2 — 6 2 ) 1 ^ 2 ,0), respectively. 
For b > s = C > the second and third terms in Eq. (0) are negative giving confocal two sheeted 
hyperboloids. Their intersections with the (rc,y)-plane and the (x, z)-plane are planar hyperbolas 
with foci at (x, y) = (±6, 0) and (x, z) = (±a, 0), respectively; they do not intersect the (y, z)-plane. 
Inverting Eq. (0) within the positive (x, y, z)-octant gives 



[X,y,Z) -\ab' bV^¥ ' av^F J {) 

with 

0<C<6<T7<a<f<l. (9) 

The remaining octants are obtained by appropriate reflections. Note that the transformation 
(x,y,z) <-> is singular on the Cartesian planes (x,y), (x,z) and (y,z), i.e. at the branch 

points of Eq. (||), see Fig. [l]. 
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Figure 1: On the planes (x,y) and (x,z) the coordinate surfaces of £, r\ and C, reduce to planar conic sections, on 
which the transformation (x,y,z) «-» (£,??, C) is singular. The figure shows them for the positive (x, y, z)-octant. 



With (p£,Pr),P<;), the momenta conjugate to (£, Hamilton's function for a freely moving 
particle in ellipsoidal coordinates reads 

- g*)(? - P) % (a 2_ ty 2 )(7y 2_ b 2 ) p g (a 2_ C 2 )(62 _ C 2 ) 3 ^ 

(£ 2 - ?7 2 )(£ 2 - C 2 ) 2 (£ 2 - ?7 2 )(r? 2 - C 2 ) 2 (£ 2 - ( 2 )(v 2 - C 2 ) 2 ' 1 ' 

The reflection at the billiard boundary £ = 1 is simply described by 

(€,V,(,P£,Pv>Pd ~* (^^C, -P£,Pr,,Pd- ( U ) 

Especially for the quantization it is useful to consider also the symmetry reduced billiard. The 
billiard is then confined to one (a;, y, z)-octant, e.g. the positive one, with the particle being 
elastically reflected when it hits the boundary or one of the planes (x, y), (x, z), or (y, z). 

The separation of Hamilton's equations in these variables can, e.g., be found in H. Because 
the Hamiltonian and the reflection condition can be separated the system is completely integrable. 
Besides the energy there are two independent conserved quantities 

K 2 = AEk = \L\ 2 + (a 2 + b 2 )p 2 x + a 2 p 2 y + b 2 p 2 z , (12) 

L 2 - ™ = *Ll + L* + #p>, (13) 
ar or y 

where L x , L y , L z denote the components of the total angular momentum L = r x p, but L ^ \L\. 
In the spherical limit a = b = we have K = \L\. Thus K is a generalization of the absolute 
value of the total angular momentum. The meaning of L becomes clear in the limiting cases of 
rotationally symmetric ellipsoids. In the oblate case (b = 0) L is the angular momentum about the 
shorter semiaxis, L = L z . In the prolate case (a = b) L is related to the angular momentum about 
the longer semiaxis, L 2 = K 2 — 2Ea 2 — L 2 . 

After separation the squared momentum can be written as 

2 _ S 4 - Iks 2 + I 

with s £ {£,,r),(}- It is convenient to take the turning points s% and S2 of ^7, C) t° parametrize 
the possible values of K and L, such that 

s 4 -2ks 2 + 1 = (s 2 - s 2 2 )(s 2 - s 2 ^. (15) 

Note that the transformation from k, I to Sj is singular for s\ = S2- 

In order to ensure real valued momenta for some configuration (£, rj, Q, Equations (|14|) and (^) 
give the conditions 

si < a 2 , < si < a, b<s 2 < 1. (16) 
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Figure 2: Bifurcation diagram (thick lines). The four regions correspond to smooth two parameter families of 3-tori, 
see Fig. ^. The meaning of the dashed lines c\ and c„ and the points Pa and P„ is explained in the text. The angle 
ip is considered in Section H. 



As for billiards in general, the energy dependence can be removed by a simple scaling, see Eq. ([14]). 

The bifurcation diagram of an integrable system shows the critical values of the energy momen- 
tum mapping from phase space to the constants of motion. Typically the critical values correspond 
to the double roots of a certain polynomial, and the different types of motion correspond to the 
ranges of regular values of the energy momentum mapping. 

In the ellipsoidal billiard the type of motion is determined by the ordering of the numbers b, 
si, a and S2- Equality in Eq. ([u]) gives the five outer lines of the bifuraction diagram, while the 
lines s\ = b and S2 = a give the inner lines, see Fig. ^. The bifurcation diagram divides the 
parameter plane into four patches. In Fig. || the corresponding types of 3-tori are represented 
by their caustics, i.e. by their envelopes in configuration space. The ellipsoidal boundary itself is 
usually not considered as a caustic. The caustics are pieces of the quadric surfaces in Eq. ([?]). 
Motion of type O is purely oscillatory in all variables £)• The oscillations in the ellipsoidal 
direction £ is given by reflections at the boundary ellipsoid £ = 1. r/ and £ oscillate between their 
caustics. The remaining types of motion are best understood by considering the two limiting cases 
of rotationally symmetric ellipsoids. Type P involves a rotation about the x-axis described by the 
coordinate rj. £ now oscillates between the caustic and the boundary ellipsoid. £ oscillates between 
its caustics. This is the only generic type of motion in prolate ellipsoids. Motion types OA and OB 
both involve rotations about the z-axis, described by the coordinate £. They are the two generic 
types of motion in oblate ellipsoids. For OA £ oscillates between the boundary ellipsoid, for OB £ 
oscillates between the caustic and the boundary ellipsoid. The way r\ oscillates between its caustics 
is different in the two cases. Motion type O can only occur in the general triaxial ellipsoid without 
any rotational symmetry. A given value of the constants of motion (E,K 2 ,L 2 ) or (E, sf,^) m 
region O corresponds to a single 3-torus in phase space. In all the other regions there exist two 
disjoint tori in phase space which have the same constants of motion. They just differ by a sense 
of rotation. The non-generic motions on lower dimensional tori corresponding to the critical lines 
in Fig. || are discussed in detail in || 10 1. 



The description of the free motion inside the ellipsoid in terms of the phase space variables 
C>P£>P»?>Pc) is rather complicated because of the change of coordinate sheets each time a 
boundary of the intervals in Eq. (||) is reached. Upon crossing one of the Cartesian coordinate 
planes (x,y) or (x,z) one of the momenta p^, p v or p^ changes from ±oo to ^oo, see Eq. (|l4|) 
and Fig. |. The singularities in Eq. (|l4|) can be removed by a canonical transformation. The new 
coordinates are better suited for the semiclassical considerations in Section || 



Triaxial Ellipsoidal Quantum Billiards 



7 




Figure 3: Caustics and boundary ellipsoid of the four types of invariant 3-tori for the ellipsoid with constants of 
the motion a) (sf , s|) = (0.05, 0.4), b) (s?, s|) = (0.25, 0.4), c) (s?, s?,) = (0.05, 0.8), d) (sj, s|) = (0.25, 0.8). 
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For the generating function of this canonical transformation we choose the ansatz 

f 2 = KOvx + Kv)ph + HOpv (17) 

The index 2 indicates that this is a generating function of type 2 in the notation of H. Goldstein, 
see p^]. Then 

, of 2 dF 2 dF 2 

x= w^ = w v= w» (18) 

are the new coordinates with (p\, p„, p u ) the conjugate momentum variables. The transformation 
is completed by relating the old and new momentum variables: 

dX d/j, dv 

Pi = ^ ft, = Tr P„ Pc, = T( P V - (19) 



(20) 
(21) 
(22) 



To remove the singularities in Eq. (|1J) we require the above derivatives to be 

dX a 



de V(£ 2 -a 2 X£ 2 -& 2 )' 

dfi a 
dv a 



dC y / {a 2_ C 2 ){b 2_ C 2 ) 



Note the negative sign of the derivative dfx/drj. These equations involve square roots of fourth order 
polynomials, i.e. they lead to elliptic integrals. Their inversion leads to elliptic functions. One finds 

«A) = „*g4 (23) 
cn(A,g) 

ri(n) = adn(^q'), (24) 

C(i>) = bsn{u,q), (25) 

where sn(</>, q), cn(<f>,q) and dn(</>, q) are Jacobi's elliptic functions with 'angle' (p and modulus q 
p7fl . Here the modulus is given by q = b/a. q' = (1 — q 2 ) 1 / 2 denotes the conjugate modulus. This 
is the standard parameterization of the elliptic coordinates by elliptic functions, see e.g. [28]. For 
the momenta one finds 

P ] = a~ s ^{s\s)-2ks 2 {s)+i), (26) 

with s G {X,fj,,v} and s(s) G T)(fJ,), ((v)} from Equations (|2^)-(l25|). The coefficients a§ are 
the signs a\ = a v = + and a^ = —. 

Transforming the coordinate ranges in Eq. @ for the old coordinates (£, rj, £) to the new coor- 
dinates gives 

0< A <F( X ,q), (27) 

0< \l <JC(q')=)C'(q), (28) 

< v < K(q), (29) 

for the motion in one octant. Here J-{Xi q) denotes Legendre's incomplete elliptic integral of first 
kind with amplitude % an d modulus q [p7], |29|. The amplitude is given by 



sin 2 X = \_ a b 2 ■ (3()) 



K,(q) is the complete elliptic integral of first kind with modulus q and K!(q) = IC(q') its complement. 
In the following we will omit the modulus in the notation for elliptic integrals because the modulus 
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will not change in the course of this paper. The appearance of the incomplete integral is due to the 
fact that we cut off the coordinate range in the ellipsoidal direction, i.e. to the billiard character of 
the underlying motion. In terms of the Cartesian coordinates the coordinate ranges in Equations 
(p7|)-(p9D yield the octant x, y, z > within the ellipsoid. Inserting (A, fi, u) into the expressions for 
the Cartesian coordinates in Eq. (||) gives 

dn(A,g) dn(/x, q') sn{v, q) 

X — a — r , [61) 

cn(X,q) 

i cn(^,g / )cn(^, q) 

V = qa — , (32) 

cn(A,<?) 

, sn(\,q)sn(u,,q')dn(v,q) 

z = qa — r . (33) 

cn(A, q) 

The functions sn(</>) and cn(</>) both have period 4/C on the real axis, dn(<j>) has period 2/C. 
Extending the ranges in Equations ( p8|) and (^9|) to the full real axis for fi and v thus gives x, y and 
z as periodic functions of fi and v. If in addition to that we let A vary in the interval [— J~(x) > J~(x)] 
the billiard dynamics becomes smooth across the planes (x,y), (x,z) and (y,z). We thus have a 
coordinate system that both separates Hamilton's equations and the reflection condition and yields 
smooth dynamics inside the ellipsoid. The motion is thus best described as a geodesic flow on the 
product of an interval and a 2-torus, 

(^,f)£[-^),%)]xT 2 , (34) 

i.e. on a solid 2-torus as depicted in Fig. ^. The flow is smooth except for the reflections at the 
boundaries A = ±J-(x) which are still desribed by the sign change 

(\V,v,V\,V^Vi>) (\V,v,-p\,p^p v ). (35) 

The whole torus 

-^(X) < A <T{ X ) , (36) 

0< n <4JC', (37) 

< v < 4/C (38) 

gives a fourfold cover of the interior of the ellipsoid. In Fig. |H| we represent the solid torus as 
a cube and mark the boundaries between the preimages of the different (x, y, z)-octants. Each 
(x, y, z)-octant gives a small cube 

[0, ±Hx)\ x [n^K!, (nf, + 1)/C'] x [n v lC, (n v + 1)/C] (39) 

with (nu,njy) G Z 2 . The fact that each of the small cubes has to be bounded by 5 neighbouring 
small cubes to make the dynamics smooth can be understood in terms of the old variables £, ry 
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Figure 5: Representation of the solid 2-torus (Fig. ^]) of the configuration space as a cube with periodic boundaries 
in the directions of /x and v. Each small cube represents one (x, y, z)-octant. They are labeled in a 'binary' way with 
respect to the signs of x,y and z, i.e. ( — , — , — ) corresponds to 0, (— , — , +) corresponds to 1, (+, + , +) corresponds 
to 7. The labels are put on the right side of each cube. 



and (. Each (x, y, z)-octant is bounded by five singular sheets of the coordinates (£, ?7, C)> see 
Fig. |l|. Note that instead of considering the three real Equations (|23|)- (|25|) it is equivalent to 
consider only the £ equation but for complex v in the fundamental domain and use the idendities 
sn(u + K. + ifc', q) = q~ l dn(u, q)/ cn(u, q) and sn(— iu + K. + iKJ ', q) = q~ l dn(«, q'). 

The four covers of the ellipsoid are related by the group of involutions which leave the Cartesian 
coordinates in Equations (|3l|)-(33) fixed. This group has three non-trivial elements 



Sx(X,fi,u) 
S 3 (X,jj l ,u) 



(-A, -\i,v) , 

(A, -n - 2K', 2K - v) 

{-\ H - 2K' ,2K - v) 



(40) 
(41) 
(42) 



Any two of them generate the group which is isomorphic to the dihedral group D 2 (also called 
"Kleinsche Vierergruppe" pOfl), 

Inspection of Equations (|31[)-(^3[) shows that it is justifiable to think of /i as a kind of rotational 
angle about the x-axis. In the y-component and ^-component /i appears as the argument of the el- 
liptic functions sn and cn which are similar to the trigonometric functions sine and cosine. Similarly 
v can be considered as a rotation angle about the z-axis. The types of motion in Fig. |3| therefore 
have the interpretations of //-rotations for type P and i^-rotations with different //-oscillations for 
types OA and OB. 

In Fig. [| each column 



-Hx),Hx)\ x [nJC', (n M + !)£'] x [0,4/C] 



(43) 



with n^eZ fixed gives a single cover of the interior of the ellipsoid. This does not hold analogously 
for v. This is familiar from the polar coordinates of the sphere where v should be compared to the 
azimutal angle and [i is similar to the polar angle. 
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For the semiclassical quantization in Section [| it is helpful to deal with a simple kinetic-plus- 
potential-energy Hamiltionian. We therefore write Eq. ( |26| ) in the form 



Ps 



+ Vg(s) 



(44) 



with 



and 



a 2 



V,(s) 



E 



(s 4 (s)-2ks 2 (s)). 



(45) 



(46) 



The effective potentials and V v are periodic functions with periods 2/C' and 2/C, respectively. 
V\ is symmetric about 0. The number of potential wells per period changes across the lines 
2k = s\ + s| = 2o 2 and 2k = s 2 + s% = 2b 2 indicated as the dashed lines c\ and c v in Fig. |2[ 
Between c\ and c u , has two maxima per period at integer multiples of K' , V u has one maximum 
per period at odd integer multiples of JC and V\ has a single maximum at A = 0. The effective 
potentials and energies for this region in Fig. ^ are shown in Fig. ||. Above c\, V\ has a minimum 
at A = and two symmetric maxima, and has only one maximum per period at odd integer 
multiples of K! . It is easy to check that the effective energy E\ is always less than the potential 
energy at the minimum at A = 0. This minimum thus has no consequences for the classical 
dynamics. At c\ in Fig. V\ changes from one to two maxima, and at P\ we additionally have 
V\(0) = E\. E^ reaches its minimum value relative to here. Below the line c v the maxima 
of V u at odd integer multiples of K, change into local minima and V v has two maxima per period. 
The maxima of at odd integer multiples of KJ have vanished here and has one maximum per 
period. Again it is easy to check that the effective energy E v is always less then local minima of V u 
at odd integer multiples of /C. The local minima thus do not influence the classical dynamics. At c u 
in Fig. ^, V u changes from one to two minima and at P u we additionally have V u ((2n + 1)/C) = E v 
for n G Z. E^ reaches its minimum value relative to here. These cases are summarized in Fig. 0. 



3 The Action Integrals 

For the calculation of actions it is useful to inspect the caustics in Fig. ||. The action integrals are 
written in the form 

1 / , m s 



Ps ds 



(47) 



with (s, s) £ {(A, £), (/i, r/), (y, £)}. The integers m s and the integration boundaries s_ and s + 
can be found in Tab. ffl, see S], also the final comments in Section For the symmetry reduced 



type 




rrir, 




£- 


e+ 


V- 


V+ 


C- 


c+ 





4 


4 


4 


a 


1 


b 







Sl 


OA 


4 


2 


±4 


a 


1 


Sl 


S2 





b 


P 


2 


±4 


4 


S2 


1 


b 


a 





Sl 


OB 


2 
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±4 
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Sl 
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Table I: Integration boundaries s_ and s+ and multipliers m s in Eq. (^) for the four types of motion O, OA, P 
and OB. 



ellipsoidal billiard any motion is of oscillatory type always giving m s = 2. To distinguish the 
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Figure 6: Effective potentials Vs (solid) together with the effective energies Eg (dotted) for the parts of the regions 
O, OA, P and OB with 2b 2 < s\ + s% < 2a 2 . At the top the reflections are indicated which determine the symmetry 
of the wave functions according to the parity at the reflection point, see Section m 
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Figure 7: Analogue of Fig. |^ for parameter combinations in Fig. ^| above the line c\ (first row), at Pa (second row) 
below the line c„ (third row) and at P v (last row). 
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Figure 8: Energy surface and separatrix surfaces in the space of the actions J. 



symmetry reduced actions from the actions of the full ellipsoid we write the former with tildes, i.e. 

1 r s + 

h = Is = - / Ps ds. (48) 

7T Js_ 

The presentation of the energy surface H{I) = E in the space of the actions / is not smooth 
because an action variable can change discontinuously upon traversing a separatrix. In contrast to 
that the symmetry reduced system H (J) is continuous. For the quantum mechanical considerations 
it is advantagous to have a continuous energy surface even for the full system. We therefore 
introduce the actions 

J = 21, (49) 

which have the property that the phase space volume below the energy surface H(J/2) = E in the 
space of the actions J is equal to the phase space volume below the energy surface H{I) = E in 
the space of the actions i" for the same energy E. In Fig. || H(J/2) = E is shown together with 
the separatrix surfaces sf = b 2 and s 2 = a 2 . Because the action variables scale with the energy the 
separatrix surfaces are foliated by rays through the origin. They divide the action space into the 
four regions corresponding to the different types of motion O, P, OA and OB. 
Inserting the momenta from Eq. (^) and substituting z = s 2 in Eq. (|47| ) shows that the action 
integrals are of the form 

r\z-s 2 )(z-sl)^- (50) 

with 

5 

w 2 = l\(z - Zi ) (51) 
i=i 

where z~, z + are sucessive members of {zi, Z2, Z3, 24, z$, Zb} = {0, s 2 , b 2 , s?,, a 2 , 1}. z\, = 1 corre- 
sponds to the boundary of the billiard. The differential dz/w has the six critical points Zi, ..-,25 
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Figure 9: Riemann sphere C with the critical points zi, ze- The rectangular boundary with the point z% should 
be considered as the point oo. The dashed line marks the boundary ellipsoid extented to the complex plane. The 
points —i and i are marked for reasons of orientation. 




Triaxial Ellipsoidal Quantum Billiards 



15 



and zq = co which implies that the integrals in Eq. (^) are hyperelliptic. There do not exist tabu- 
lated standard forms for these integrals but there is the well developed theory of so called Abelian 
integrals. The main object of this theory is a Riemann surface, in our case the hyperelliptic curve 

5 

TZ w = {(z,w) eC 2 :w 2 = H(z- Zi )}. (52) 

i=l 

Here C denotes the compactified complex plane, i.e. the Riemann sphere. To construct a picture 
of 1Z W we proceed in the following manner. We order the critical points Zi, i = 1,...,5, according 
to their magnitudes. This gives 4 orderings, one for each type of motion, e.g. Z\ = 0, Z2 = s\, 
Z3 = b 2 , Zi = s 2 , and z§ = a 2 for type O, see Tab. |. The points z% are marked on the Riemann 
sphere, see Fig. |9[ We then slit the Riemann sphere along the real axis between the points Z{ and 
Zi+i for i = 1,3, 5. Excluding the three slits from C the sign of w is everywhere well defined on this 
manifold when it is fixed at one arbitrary point. In Fig. ^ we choose the sign of w to be positive 
right above the slit [z\, z<i\. Then the sign is negative right above the slit [z^, 24] and again positive 
right above the slit [z&, zq]. Right below the slits the sign of w is opposite to the sign right above. 
Around the slits we have the closed paths z\ ^± Z2 , Z3 ^ z<± and z$ ^ zq. On another copy of C we 
introduce the same slits but choose the sign of w opposite to the choice on the former copy. The 
path from Z2 to Z3 on the former copy in Fig. || is assumed to be the first half of a closed path 
Z2 ^ Z3 of which the second half from Z3 back to Z2 lies on the latter copy. The same is assumed 
to hold for the closed path 24 ^ Z5. To unify the view glue the two copies at the corresponding 
slits such that the corresponding critical points coincide and such that w changes smoothly across 



the seams, see Fig. |10|. The result is a compact Riemann surface, i.e. a manifold which carries a 
complex structure and to which the full machinery of Cauchy integration theory is applicable. The 
surface has genus g = 2 and there are 4 non-contractable paths on the manifold which cannot be 
transformed smoothly into each other. They form a basis of the four-dimensional homology group 
corresponding to this surface. The homology basis may be specified by the choice of the closed 
paths zi ^ Z2, Z2 ^ Z3, Z3 ^ Z4 and 214 ^ 25. The path z§ ^ zq is homologous to the sum of 
z\ f± Z2 and Z3 ^ Z4 then. From the non-trivial topology of the Riemann surface 7Z W it follows that 
there may exist non- vanishing closed integrals (even for vanishing residues). The action integral 



Eq. (50) is of this type. It is an integral with singularities but vanishing residues - a so called 
Abelian integral of second kind. The actions integrals I u and 1^ of Eq. (|50| ) are taken along the 
closed paths z\ ^ Z2 and Z3 ^ z^. Due to the reflection at the boundary ellipsoid the action 
integral I\ is not taken along a closed path. It is taken along the slit [z$, zq], but only between z$ 
and Zf,. It is therefore called incomplete. The integrals 1^ and I u are called complete. These three 
integrals give real numbers. In contrast to this the integration of Eq. (^) along the closed paths 
Z2 ^ Z3 and Z4 f± 25 yields purely imaginary numbers. These integrals have an important physical 
meaning for the semiclassical quantization scheme in Section ||. They give the penetration integrals 
which will be needed for the discussion of quantum mechanical tunneling. At this stage we already 
mention that there are only two such penetration integrals and we define them as follows: 

@u = e c = -2i I p c d(, (53) 

J Sl 

@x = e ( = -2i Hptd^. (54) 

J a 

The factor i in the definition turns both integrals into real numbers. It is useful to take these 
definitions independent from the type of motion O, OA, P and OB, i.e. for any ordering of b, s\, 
a and S2- This will become clear in Section ||. 



In Fig. we show the ranges for the coordinates s 2 G {£ 2 , rr 2 , £ 2 } on the circle in the parameter 
plane in Fig. ||. Generically the ranges for £ 2 and rj 2 and the ranges for r] 2 and £ 2 are separated 
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Figure 11: The parameter ranges for s 2 £ {C 2 ! 7 ? 2 ^ 2 } define a family of hyperelliptic curves parametrized by the 
angle ip in Fig ^| 



by finite gaps. This gives the Riemann surfaces of genus 2 as described above. Now consider the 
circle in Fig. 0. On the bifurcation lines which are reached for angle ip equal to integer multiples 
of 7r/2 one of the two gaps vanishes. One of the penetration integrals in Eq. ( |53| ) then vanishes 
too. This means that two of the three slits in Fig. |9| merge and the genus of the Riemann surface is 
diminished by one. On the bifurcation lines we thus find elliptic curves (genus 1) with an additional 
pole in the differential for the actions. They can also be considered as singular hyperelliptic curves. 
The three action integrals and the one remaining penetration integral are of elliptic type then. In 
[||] analytic expression in terms of Legendre's standard integrals are calculated for these cases. 

4 The Quantum System 

The quantum mechanical billiard problem is the problem of determining the spectrum of the Lapla- 
cian in the billiard domain with Dirichlet boundary conditions on the boundary. Equivalently, this 
is given by Schrodinger's equation for a free particle in the ellipsoid which in turn is Helmholtz's 
equation in three dimensions, 

- = Eip. (55) 

As in the classical case the potential vanishes inside the ellipsoid and is infinite outside the ellipsoid. 
This potential classically leads to elastic reflections and quantum mechanically imposes Dirichlet 
boundary conditions on the ellipsoid. The three discrete symmetries of the ellipsoid are the reflec- 
tions at the three Cartesian coordinate planes. The wave function can have even or odd parity with 
respect to each discrete symmetry, ip(x,y,z) = TT x ip(—x,y,z) etc. Combining the two parities for 
each dimension we obtain a total of eight parity combinations denoted by tv = {^x^y^z) where 
each parity is from {+,—}• 

Corresponding to the two sets of classical coordinates we get two sets of quantum mechanical 
equations. In both cases the separation is the same as in the classical case and the wave function 
ip is a product of three separated wave functions. The ellipsoidal coordinates (£, rj, £) lead to the 
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analogue of Eq. (|l4| ) , which is 

-h 2 (V( s2 -a 2 )(s 2 ~b 2 )-^j Ms) = 2E (s 4 - 2ks 2 + l) i/j s (s) (56) 

with s £ {£, 77, £}. If we set E" = but keep K and L finite we obtain one of the many forms of the 
Lame equation (2^, |3l] ]. Since we are not only interested in the solution of the Laplace equation 
in the ellipsoid but in the spectrum of the Laplacian we have to consider this generalized Lame 
equation, known as the ellipsoidal wave equation. 

Transforming the equation into the regularized coordinates leads to the analogue of Eq. (p6|) . 



if 

- h 2 a 2 a s —^ s (s) = 2E (s(s) 4 - 2ks{sf + I) ^ s (s) (57) 

where s £ {A,/j, v} and s(s) G {^(A), 77(/x), C( t/ )} from Equations (p3|)-(p5|). Comparing to the 
equations for the billiard in the ellipse |32| Eq. ([5?]) is analogous to the Mathieu equation(s) in its 
standard form, while Eq. (56) is analogous to its algebraic form. Note that similar to the classical 



case it would be sufficient to only consider the equation for v in the complex domain instead of the 
three equations for real arguments. 

The Dirichlet boundary conditions require that the wave function ?/>(A, /U, u) = vp\(\)'ip^(n)'ij; l ,(i>) 
is zero on the ellipsoid, which gives ip\[±J-(x)) = 0. The solutions in the two angular variables 
and v must be periodic with periods 4/C' and 4/C, respectively, in order to give a smooth function 
on the solid 2-torus described in Section ^. 

For \x and z/, Eq. (|57|) is a linear differential equation with periodic coefficients. Floquet theory 
guarantees the existence of solutions tp^ with period an integer multiple of 2/C' and solutions ip v 
with period an integer multiple of 2/C, respectively. The involutions Si, S2 and S3 in Equations 
(|40|)-(p2|) relate the symmetries of the separated wave functions to the parities ir x , 7r y and tt z . 
Starting with S± the separation of the invariance condition ^(Aj/i, v) = ip(Si(\, fj,,v)) gives 

V>a(A) _ ^ m (-m) ( 58 ) 



Since the left hand side and the right hand side are functions of A and /x alone they have to be 
equal to some common constant. Because we may change the sign of A and /i independently giving 
the reciprocals of both sides of Eq. (|5^ ) the separation constant must have unit modulus. From 
Eq. ( |33| ) we see that the sign is the parity ir z . Similarly, from the invariance of tp under S2 we get 

Vv(-/.-2/C') = ijjjy) ^ 

From replacing [i by —\i — 2K! and/or v by 2/C — v we see that both sides of Eq. (|59|) again have to 
be equal to a separation constant of unit modulus. With the aid of Eq. (32) we may identify the 



sign with the parity n y . ir y thus gives the parity of the wave function tp^ for reflections about /C' 
and of ip u for reflections about /C. From the invariance of ip with respect to S3 we get 

V> A (-A) _ M") (m) 

V^a(A) ^ v {2K,-v) 1 > 

or with the results from above 

^(/J - 2/C') = ITyTTzlp^fJ,) . (61) 

This relates the product of the parities ix y and tt z to the period of tp^. tpn is 2/C-periodic for 
%7Tz = + and 4/C'-periodic (i.e. not 2/C-periodic) for TT y ir z = — . Similar arguments hold for the 
wave function tp u . Here tt x gives the symmetry of ip v for reflections about 0. The product of the 
parities tt x and ir y determines its period: tp v is 2/C-periodic for n x ii y = + and 4/C-periodic (i.e. not 
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2/C-periodic) for 7r x TT y = —. The parities for the separated wave functions are shown at the top of 
Fig.g 

Even though the ellipsoidal coordinates (£, 77, Q are not regular, the parities are most simply 
expressed by properties of the wave functions in these singular coordinates. Let us first rewrite 
Eq. ( |56|) in the form 

ftp" + gtp' + htfj/h 2 = (62) 



with the polynomials 



f(s) = ( s 2 -a 2 )(s 2 -b 2 ), (63) 
g(s) = f'( s )/2 = s(2s 2 -a 2 -b 2 ), (64) 
h(s) = 2E{s A - 2ks 2 + I) . (65) 



The singularities of Eq. (p2\ j and equivalently of Eq. flog ) are given by the zeroes ±a and ±6 of 
/. We postpone the question of additional singularities at infinity to Section |6| because they are 
not important for our numerical calculations. In order to look at the asymptotics of the solutions 
of Eq. ( |6"2| ) at the singular points we calculate the corresponding indicial equations. Denoting 
the position of the singularity under consideration by c, the exponents a of the solutions are the 
solutions of the indicial equation (see e.g. |28|]) 



a 2 + ( Pc - l)a + q c = 0, (66) 

where 

Pc = lim 77~|( s ~ c )' ?c = lim-^^(s-c) 2 . (67) 

To calculate p c it is best to perform the partial fraction decomposition of g/ f, 

g(s) 1 / 1 1 1 1 \ 

' + 7 + + • (68) 



f(s) 2 \s — b s + b s — a s + a 

From this it is obvious that p c = 1/2 for all singular points. Since h/f only has simple poles q c = 0, 
the exponents are and 1/2. Because the singularities of the wave equations do not produce 
essential singularities in its solutions these singularities are called regular. In our case the two 
exponents refer to the two parities possible at a regular singular point. We require ip v (b) = V'c(^) = 
for TT y = — and ipnQ>) = V'c(^) = 1 ( U P t° normalization) for ir y = +. Similarly the value at the 
regular singular point a of the wave functions ip^ and ip v determines the parity ir z . For ir x it is 
a little simpler, because it is determined by the value of a t the ordinary point ( = 0. The 
boundary condition at £ = 1 always is V^(l) = 0. The need for the solution to be invariant under 
an additional symmetry group (arising e.g. if we work on a covering space) does not appear, because 
we only solve the wave function in one octant. The boundary conditions are summarized in Tab. |ll[ 
Note that in line with the above considerations the table shows a very simple structure: the sign — 
or + in the first three parity columns successively determine the value or 1 of the wave functions 
at 0, b and a. 

In our numerical procedure we are going to start integrating at the regular singular points. 
Since this is impossible for initial conditions belonging to the solution with exponent a = 1/2 we 
have to factor out this behaviour analytically. To find solutions with the parities ir z = — and/or 



we employ the transformations 



TT y = +,7T 2 = - 
TTy = — ,7T 2 = + 
7T„ = — , 7t z = — 



tjj = \J s 2 — a 2 ip^ , (69) 

ijj = Vs 2 - & 2 V>-+, (70) 
if, = ^s 2 -a 2 Vs 2 -b 2 ?p—, (71) 
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Table II: Parities, boundary conditions of the separated wave functions ip B , periods of the separated wave functions 
ip v and and Cartesian nodal planes. 



and leave ip = unchanged for TT y = ir z = +. The polynomials h and g in Eq. (^) change 

according to 

h + _(s) = h(s) + 2s 2 - b 2 , g+~(s) = g(s) + 2s(s 2 - b 2 ), 

h-+(s) = h(s) + 2s 2 -a 2 , g- + (s)=g(s) + 2s( y s 2 -a 2 ), (72) 

h__(s) = + 6s 2 - b 2 - a 2 , g—(s) = g(s) + 2s(2s 2 - b 2 - a 2 ). 

The functions h ++ = h and g ++ = g remain unchanged. The resulting transformed equations 
change the prefactor 1/2 in Eq. (|68| ) to 3/2 for the terms involving ±a, ±6, or both, respectively. 
Hence p c = 3/2 at the correponding regular singular point and a = 0, —1/2. We are now able to 
start integrating at the singular points c = a or c = b, or to be more precise, a distance As away 
from them, always with the special velocity that corresponds to the regular solution with exponent 
a = 0. The initial conditions are 

t//(c± As) = -^4, $(c± As) = 1± As^'(c± As). (73) 
5(c) 

In order to find t/j, three conditions on the three separated wave functions ip s have to be fulfilled 
simultaneously. This is possible because there are three parameters E, k and I in the three equa- 
tions. However, each equation depends on all the three separation constants; the equations are 
separated but the constants are not. We use a numerical procedure similar to that described in 
p2| , the essential difference being that for the ellipsoid the wave function ip^ has a regular singular 
point on both ends of the interval. Since it is not possible to integrate a regular solution into a 
singular point, but only away from it, we divide the interval into two equal parts [b, (a + b)/2] and 
[(a + b)/2, a] and require the solution to match smoothly at s = (a + b)/2. This is called shooting to 



a fitting point [33]. These two and the two remaining intervals [0, b] and [a, 1] are all transformed to 
[0, 1], and the resulting system of four equations is simultaneously solved. With Newton's method 
|33|| the three free parameters are adjusted to satisfy the three remaining conditions = 0, the 

smoothness condition at the fitting point (a + b)/2 for ip v and V^(0) = f° r 11 x = + or ^(0) = 
for tt x = — , respectively. Taking the semiclassical values for E, k and I from Section || as an initial 
guess, the method always converges to the exact eigenvalues. 

Because we are free in the normalization of the three separated wave functions ip s they can 



be multiplied by constant factors to give one smooth function on the interval [0,1], see Fig. 12 
ip has n zeroes 6 (0,6), m zeroes E (6, a) and r zeroes S (a, 1). The quantum numbers (r,m,n) 
together with the parities it = (TT x ,7r y ,ir z ) completely determine the state which we denote by 
\r,m,n;ir x Tr y ir z ) . The quantum numbers (r,m,n) belong to the reduced system in one octant; the 
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Figure 12: The eight transformed eigenfunctions tjj with all quantum numbers 1 and all possible parities, i.e. 
|1,1,1;±±±). The parity n x is — for wave functions starting at and + otherwise. The remaining two parities 
ir y ,n z are distinguished by the dashing. 



parities determine the wave function on the boundaries of the octant. It is not so simple to count 
the corresponding number of nodal surfaces in the full system. It is complicated by the fact that 
the (x,y)-plane and (x,z)-plane are composed of two different types of quadrics, see Fig. |l|. Away 
from the three Cartesian coordinate planes the number of ellipsoidal nodal surfaces is counted by 
r, the number of one sheeted hyperboloidal nodal surfaces ("rotating" about the shortest semiaxis 
z) is given by m and the number of two sheeted hyperboloidal nodal surfaces ("rotating" about 
the longest semiaxis x) is given by n, but because each surface has two sheets this makes In nodal 
surfaces. Depending on the parity combinations tv the Cartesian coordinate planes give additional 
nodal planes according to the last column of Tab. ||. 



5 Semiclassical Quantization 

The semiclassical quantization of the ellipsoidal billiard is obtained from single valuedness condi- 
tions that are imposed on WKB wave functions on the fourfold cover discussed in Section Let 
us consider Schrodinger's equation for a general one dimensional Hamiltonian of the form 



H=--^ + V(q). (74) 



For a fixed energy E in each region j between two successive classical turning points a WKB wave 
function of the form 

V { %) = (4° eMiS,(q)/h) + A® eM-iSM/V) / VpW) (75) 

is reasonable. Its phase is given by 

Sj(q) = [ 9 p(q')dq' (76) 



with the classical momentum p(q) = \Jl{E — V{q)). A^} and A_ are constants. The reference 
point <7j for the phase integral is an arbitrary point in the region under consideration but it is 
convenient to take it as the left or right classical turning point although the WKB wave function 
is a good approximation only away from the classical turning points. 
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Figure 13: Effective potential (solid line) and effective energy (dashed line) for motion type OA. The turning 
points fii (i — 1, 8) define WKB wave functions V'm ' whose amplitudes are connected by the matrices P, Mi and 
M 2 . 



In the following we will consider WKB wave functions only in classically allowed regions although 
they are valid even in regions where E < V(q) giving real exponentials in Eq. (^). The amplitude 
vectors A (1) = (A^jA^'f and A (2) = {A^^A^f of WKB wave functions in two classically 
allowed regions 1 and 2 separated by a classically forbidden region are related by the matrix 
equation = M(@)A^ with the tunnel matrix (see |34| and the references therein) 



M(0) 



e/h 



y/l + e -2S/h 



y/l + e -2B/h 



(77) 



where 

= -i p{q) dq (78) 
Jqi 

is the penetration integral of the potential barrier. Here q\ and q2 are the turning points to the 
left and right of the barrier, i.e. V{q) < E for q < q\ and q > q2 arid V(q) > E for q\ < q < q^. 
The matrix (|7?]) remains valid if we increase the energy E above the barrier's maximum. Then the 
classical turning points become complex (gi complex conjugate to (72) giving a negative penetration 
integral in Eq. (f7^). For —0 3> h the matrix M(O) becomes the identity matrix. 

The amplitude vectors A^ and A^ of two WKB wave functions defined in the same classically 
allowed region but with different reference points for the phase integral are related by the phase 
shift A^ = P{4>)A^> with the matrix 



exp(i(j)/h) 

exp(—i(f)/h) 



(79) 



where <j) = Jq^ p(q) dq. 



Let us now specify the Hamiltonian ([74]) for the ellipsoidal billiard by the consideration of 
the effective potentials and energies defined in Equations (|44|)-(|4"6|). To illustrate the semiclassical 
quantization scheme we concentrate on the [i degree of freedom and again present the effective 



potential and energy E^ for motion type OA in Fig. 12 . In the range [0, 4/C'] we have the eight 
turning points /Xj marked in the figure. Taking them as the reference points for the definition of 
WKB wave functions we get two wave functions in each of the four classically allowed regions. The 
amplitude vectors A^ and A^ are connected by the phase shift matrix (|79|) with <j) = p^ dfi. 
From Tab. | and the negative sign in Eq. (|2l| ) it becomes clear that (j) - 
action defined in Eq. (|49|). The matrix 



(7r/2)J M with J M the 



P 



exp(-i(ir/2)J lx /H) 








exp(i(7r/2)J M /^) 



(80) 
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then also relates the pairs of amplitude vectors A^ and A^ , A^ and A^ , A^ and A^ , see 
Fig. 0. 

The amplitude vectors A^ and A^ are related by the tunnel matrix ( |77|) where the penetration 
integral G = — i p^ dp is the penetration integral @ v defined in Eq. (|5S|) (again see Tab. | and 
keep in mind the negative sign in Eq. (f2l|)). The integration boundaries in the definition of Q u 
were independent of the classical type of motion. Therefore the connection relation remains valid 
if the effective energy and potential change such that we classically have a different type of motion, 
especially for motion types O and P where the turning points p 2 and ^3 become complex, see Fig ||| 
We set 

Mi = M(-B v ) . (81) 

The matrix M\ also connects the amplitude vectors A^ and A^ 7 \ 

Similarly one finds that the pairs of amplitude vectors A^\ A^ and A^ 8 \ A^ are related by 
the tunnel matrix 

M 2 = M(—®\) (82) 
with 0\ from Eq. (^) where we have taken into account the 4/C'-periodicity of the wave function 

Starting at p\ the quantization of the p degree of freedom now reduces to finding an effective 
energy and an effective potential for which there exists a non-zero amplitude A^ which is 
mapped onto itself upon one traversal through the interval [0,4/C'], see Fig. This is equivalent 
to the quantization condition 

det((M 2 PMiP) 2 - 1) = (83) 



with 1 the identity matrix. Similar quantization conditions can be found in [35, 37, 38] . Because 
of detM(Q) = detP(0) = 1, Eq. Q may be rewritten as 

tr(M 2 PM 1 P) 2 = 2. (84) 

The eigenvalues of Eq. ( |83"D include all parity combinations ir y and tt z . To distinguish between the 
different parities more information is needed. The parities give the additional conditions 




TTy 



1 

1 



A<® , (85) 



These conditions have to be solved consistently with the above tunnel relations. From the various 
possibilities to do this we choose the following. We map the amplitude vector A^ from a point 
p G (pi,p 2 ) to the point p. + 2K! G (ps,pq). From Eq. (|6l|) we know that this produces the 
sign TT y TT z . There are two possibilities to replace one of the tunnel matrices in this map by the 
corresponding condition in Equations (|85|) and (jB6|) . We thus get the equations 



with the matrices 



B 



C = M 2 P7T y l j P-TT y 7T z l. (89) 



BAW=0, CA^=0 (87) 
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Eq. (|87|) is the analogue of Eq. (|83| ) for half the interval [0, 4/C']. We are free in the normalization 
of the WKB wave function. Setting = 1 we find A { y = -Bu/B 12 . If we insert this into 
the equation involving the matrix C and decompose the resulting equations into their real and 
imaginary parts the remaining independent conditions are 

(0 x +0„)/h _ i 

cos^/ft) = . z yC (90) 



^(l + e 20 A /ft) ( 1 + e 26„/ft) 



and 



sin(7r Vft) = ^ +%e . (91) 

J(l + e 2e A /R) (i + e 2e,/n) 

These equations have to be fulfilled simultaneously. They are not independent of each another, but 
the relation is simple: the second equation is fulfilled on every second solution of the first equation. 

For the A and v degree of freedom we have to comment on the additional potential barriers 
in Fig. |7] appearing below the line c v and above the line c\ of the bifurcation diagram in Fig. |2[ 
As we have mentioned in Section |2| the effective energies always lie below the additional local 
potential minima. They only reach the local minima at the points P v and P\, respectively. The 
corresponding classical turning points are always complex and therefore these barriers would enter 
the quantization scheme almost always with large negative penetration integrals, i.e. with tunnel 
matrices close to the identity matrix. The only exceptions occur in the regions close to the points 
P v and P\ which lie at the border of the bifurcation diagram. Here the action goes to zero, i.e. 
the semiclassical approximation is expected to give poor results anyway. The additional barriers 
will therefore not be taken into account. The quantization conditions for A and v are then exactly 
the same as in the case of the planar elliptic billiard discussed in |^2|. We only state the results. 
For the A degree of freedom one finds the two conditions 

cos(vr J x /h) = . (92) 



and 



-1 

y/l + e -2B x /h 



sin(vr J x /h) = - 1 (93) 



and for the v degree of freedom the conditions 



and 



cos(TrJ u /h) = - x y (94) 

sm{TrJ v /h) = === • (95) 

The actions J x and J v are again taken from Eq. ( [49"1) and the only penetration integrals that appear 
are those defined in Equations (p3) and (|54[). 



We first inspect the quantization conditions in Equations (p^)-(pB|) for the limiting cases 
I ©A |) \ ®v\ ^ ^- The signs of the penetration integrals Q x and @ u determine the type of classi- 
cal motion. The limiting cases \@ x \, 3> h thus correspond to the four regions in classical action 
space far away from the separatrix surfaces in Fig. ||. From the limiting quantization conditions for 
the actions J the quantization of the original action variables / can be deduced from Tab. |. We 
summarize the results in Tab. |III| . The limiting quantization conditions for i" may be compared with 
the EBK quantization in Eq. (||). From the identification of the Maslov phases a = (a x ,a^,a u ) in 
the equations in Tab. |TfJ we find 

ol = (4, 2, 2) for type O, a= (4, 2, 0) for type OA, (96) 
a = (3, 0, 2) for type P, a = (3, 2, 0) for type OB. (97) 
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Eq. 



Eq 



type O: -@x,@„ > K 



© 
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= Jx = 


(n\ 






I h 


= J\ = 


(»A + |)« 




o J 


™A 


= 2r + 


(1- 


t»)/2 





n\ 


= 2r + 


(l-7T,)/2 
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1 




= J M = 


(% 




-1 


1 ^ 


= -V 2 


= + 




% J 




= 2m- 




" 7T„)/2 







= m 
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1 




= Jv = 










= ±J„ 




© 






= 2n + 


(1- 


7r*)/2 







= 2n + 


(2 - 7r B - 7r y )/2 



type P: 6A,e„ > ft 



© 


°1 ' 




= Jx/1 


= {nx + \)K 







= Ja/2 


= (nx + \)h 


© 






= r 




-1 


ha 


= r 




© 




\ >» 




= ±n M fi 





1 h 




(n p + 


© 


o J 




= 2m 4 


- (1 - 7r„7r,)/2 






= 2m 4 


-(l-7T,)/2 


© 


1 




= Jv = 


(n„ + ^)7i 






= 


= ±n v h 


© 




n v 


= 2n + 


(1 - 7r x )/2 







= 2n + 


(2-n x - ir y )/2 



type OA -6a, -6^ > ?i 



type OB: O x ,-e„ > ^ 



Table III: Limiting quantization conditions for the 4 types of classical motion. The left hand sides of the braces in 
each box gives the limiting value of the right hand side of the equation cited in the very first column. The non-negative 
integers n\, n M and n v are related to the quantum numbers r, m and n introduced in Section ^, see Tab. Q. 



The Maslov indices a M and ol v are in agreement with the simple EBK rule stated in the introduction: 
For motion of type O [i and v oscillate, for motion of type P the motion is rotational in \x and 
oscillatory in v, OA and OB are oscillatory in the \x degree of freedom and rotational in see 
Section ||. For the A degree of freedom we have to take into account the reflection at the boundary 
ellipsoid which wave mechanically leads to Dirichlet boundary condition. For motion types O and 
OA A oscillates with two reflections giving a\ = 4. For motion types P and OB A oscillates between 
the boundary ellipsoid and the caustic giving a\ = 3. 

The EBK quantization condition in Eq. (|j) defines a lattice in classical action space. The 
Maslov indices determine how this lattice is shifted relative to the simple lattice (nh). Since we 
have four different vectors of Maslov indices for the ellipsoidal billiard we have four different lattice 



types away from the separatrix surfaces in Fig. [8|. We present the different lattices in Fig. 14 
for quantum cells of width AJ\ = AJ^ = AJ U = 2h. Each cell contains eight quantum states. 
For motion type O all states are non-degenerate. For motion types OA and OB each states is 
twofold quasidegenerate according to the two senses of rotation in v. Analogously for motion type 
P each state is twofold quasidegenerate according to the two senses of rotation in the variable /x. 
From the quantum mechanical point of view the quasidegeneracy can be understood in terms of the 
effective energies and potentials in Fig. ||. For eigenvalues which classically correspond to rotational 
motions far away from the classical separatrices the effective energy is much larger then the effective 
potential. The energy is then dominated by the kinetic energy, the specific shape of the potential 
becomes irrelevant. The effective energy then only depends on the net number of nodes of the wave 
function and not on the location of the nodes. Therefore wave functions with different symmetries 
but the same net number of nodes give the same effective energy. With the aid of Tab. Ill we can 
identify the states corresponding to the capital letters in Fig. 113, see Tab. |TV|. 



The quantization conditions in Equations (|92D-(|95[) are uniform, i.e. they do not only give 
the limiting EBK lattices in Fig. [14| but also specify how these lattices join smoothly across the 
separatrix surfaces of Fig. |8[ In the following we will refer to the uniform lattice in action space as 
WKB lattice. The transitions may be described in terms of effective Maslov phases. In order to 
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O OA 




Figure 14: Quantum cells A J\ — AJ F = A J„ = 2 in classical action space for the 4 limiting cases of classical types 
of motion. J is measured in units of h. 
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Table IV: Quantum states in Fig. [l4| 
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see this we insert the EBK-like quantization conditions for the symmetry reduced ellipsoid 

I = J/2 = {ii + a./4)h (98) 



with fi = (r,m,n) into the left hand sides of Equations (|92j)-(|95|). The parities n x , 7r y and tt z 
on the right hand sides determine whether we have Dirichlet or Neumann boundary conditions 
on the corresponding piece of the Cartesian coordinate plane bounding the symmetry reduced 
ellipsoid. The different parity combinations altogether give the quantum states of the full ellipsoidal 
billiard. We now solve Equations (|92"|)-(p5|) for the effective Maslov phases a. The quantum numbers 
h = (r, m, n) drop out because of the 2-7r-periodicity of the trigonometric functions and it remains 
to invert the tangent on the correct branch which is determined by the parity combination. A little 
combinatorics gives 

2 

a x = vr z -arctan e 0A//ft + 3 -vr z , (99) 

7T 

2 

a v = 7Ty — arctan g®"/^ -)- 2 — tt x — TTy. (100) 

For this simple form cannot be achieved. Instead we write 

a„ = ^arg (n z 7T y e ( - @ " +e ^ h -l + i (ir z e @ ^ h + % e e "/ a ) ) , (101) 

where arg maps the polar angle of a complex number to the interval [0, 2ir). Essentially the effective 
Maslov phases consist of the simple switching function (2/tt) arctan e x which changes from to 1 
when x changes from — oo to +oo. From the simple form of the effective Maslov phases it follows 
that the ranges in which the semiclassically quantized action variables J may vary are restricted 
according to J = (J\, J^, J u ) mod 2H £ P-tj-, where the parity boxes P^ have side length h/2 in the 
directions of J\ and J v and side length h in the direction . For the different parity combinations 7r 
we find 

P___ = [(3/2)ft, 2h] x [h,2H] x [(3/2)h,2h] , 

= [h,(3/2)H] x [(1/2)^,(3/2)^] x [(3/2)h,2h], 
P_+_ = [{3/2)h,2h} x [(1/2)^,(3/2)^] x [^,(3/2)^], 
= [h,(3/2)h] x [0,h] x [H,(3/2)h], 

= {(3/2)h,2h\ x [(h,2h] x [(1/2)M], 
P+_+ = [H,(3/2)h] x [(1/2)^,(3/2)^] x [(l/2)h,h], 
P++_ = [(3/2)H,2h] x [(l/2)h, (3/2)ft] x [0,(1/2)^], 
P+++ = [h,(3/2)h] x [0,h] x [0,(l/2)ft], 



see Fig. 15 



Let us comment on the numerical procedure to solve the quantization condition (|98|) . For 
given quantum numbers n and parity combination n we have to find the corresponding zero of the 
function 

G(E, si s 2 2 ;ir, h) = I(E, s\,sl) - (n + a(E, a?, s 2 2 ,iz)/4)h . (103) 

As in the exact quantum mechanical problem this problem is not separable for the separation 
constants and we have to apply Newton's method in three dimensions. In order to find good 
starting values for Newton's method we introduce an approximate function G app for G. The 
functional dependence of G app , i.e. of the approximate actions i" app and Maslov phases a app , 
on the parameters (E, sf,s|) should be very simple such that the zeroes of G app can be found 
analytically. For Q app we simply take the mean value of a. for a given parity combination n. In 
order to get an expression for / app we take advantage of two properties of the energy surface in 
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Figure 15: The parity boxes P-k that contain the semiclassical states. P-++ is obscured by the other parity boxes. 
J is measured in units of h. 



action space. Firstly the shape of the energy surface is very similar to a triangle and secondly up 
to a simple scaling the shape does not change with the energy. If we denote the intersections of 
the energy surface E = 1/2 with the coordinate axes in action space by I\, 1^, and I v the action 
variables can be approximated by 



I&pp(E, s±, s 2 ) 



'2E 





V / 



+ 7i 



f-h\ 




+ 72 



( -h\ 



\ /J 



(104) 



where 71 and 72 parametrize the approximate triangular energy surface E = 1/2. This can be 
considered as a crude periodic orbit quantization: we take the actions of the three stable isoltated 
periodic orbits of the system and approximate the whole energy surface by that of a harmonic 
oscillator that would have isolated stable orbits with those actions. The analogy has to be taken 
with care because our actions scale with y/E, while those of the true harmonic oscillator are linear 
in the energy. The edges of this approximate energy surface are given by 71 = 0, 72 = 0, and 
7i + 72 = 1, respectively. In order to give s\ and S2 as functions of 71 and 72 it is useful to have a 
look at the asymptotic behavior of the actions and their approximations Q104| ) upon approaching 
the edges of the energy surface. A simple calculation shows that l v behaves quadratically in s\ for 
s\ —* 0. We set 

(105) 



a 2 7i . 



A similar consideration of the asymptotics of I\ for S2 — > 1 gives 



si 



l-(l-& 2 )(l-(7i+72)) 2/3 - 



(106) 



In our numerical calculation the starting values obtained from these approximations always were 
sufficiently good to make Newton's method converge to the right state. The quasidegeneracy of the 
states is no problem here because the degenerate states correspond to different parity combinations 
7T. In Tab. [V| the semiclassical eigenvalues are compared to the exact quantum mechanical results. 
The semiclassical energy eigenvalues are always a little too low. The same is true for k sc while l sc 
tends to be too low. We do not have a good explanation for this. 



Let us first consider the four transitions of the WKB lattice in action space across the separatrix 
surfaces of Fig. || away from the intersection line of the separatrix surfaces. Upon each crossing 
only two effective Maslov phases change appreciably. We therefore take the action component of J 
corresponding to the effective Maslov phase that stays approximately constant as being semiclassi- 
cally quantized. To do so we have to fix the quantum number belonging to this action component 
and the parities appearing in its effective Maslov phase. The quantum numbers corresponding to 
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Table V: The quantum mechanical eigenvalues (-E qm , fe qm , Z qm ) and the semiclassical eigenvalues (E BC ,k sc ,l EC ) of 
the ellipsoidal billiard for the < 24 and 1000 < £ qm < 1003. The relative error — (-£/ q m — -Esc)/-£/ q m 

is given in percent. 



the two other action components and the remaining free parities then define a family of surfaces in 
action space, which intersect the plane corresponding to the action component that already fulfills 
the semiclassical quantization condition. In Fig. ^ we represent the plane corresponding to the 
semiclassically quantized action component and the intersection lines projected onto the plane of 
the two remaining free action components. The semiclassical eigenvalues appear as the intersection 
points of the intersection lines as far as they are contained in a parity box. 

Let us first consider the transition from region O to region OA in Fig. 16 a. Upon this transition 



we always have —©a 3> h, see Tab. [IF From Eq. (|99|) we see that the effective Maslov phase ot\ 



stays approximately 3 — tt z . We semiclassically quantize J\ by fixing the quantum number r = 10 
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and the parity tt z = — . The transition of the WKB lattice takes place in the action components 
and J u . In the region corresponding to motion of type O all quantum states are non-degenerate. 
Upon the transition from region O to region OA quantum states with the same product of the 
parities tt x and TT y become quasidegenerate. For (TT x ,Tr y ) = (±,=f) the quasidegenerate states have 
the same quantum numbers (r, m,n); for (ir x ,n y ) = (+,+) and (ir x ,TT y ) = (— ,— ) they differ by 1 
in the quantum number n, see Tables [III] and [IV| and Fig. |14|. The picture for ir z = + is similar 
and therefore is omitted. 

In Fig. ^b the transition from region O to region P is presented. Here we again have a transition 



from non-degeneracy to quasidegeneracy. From Tab. [in] we see that we always have Q u h giving 
a v ~ 2 — ir x , see Eq. ( |100| ). For the semiclassical quantization of J u we choose n = 5 and tt x = — 
and we represent the transition of the WKB lattice in the components J\ and J^. In region P 
the quantum states with the same product of parities ir y and tt z are quasidegenerate, again see 



Tables III and IV and Fig. 14. For ir x = + the picture is similar and therefore is not shown here. 

For the transition from region P to region OB we always have 0a ^> h giving a\ ~ 3. The 
transition of the WKB lattice takes place in the components and J u , see Fig. [l6| c. For the 
semiclassical quantization of J\ we have chosen r = 5. Upon the transition the quasidegeneracy in 
region P explained above changes to the quasidegeneracy in region OB. Here quantum states with 
the same product of the parities ir x and ir y are quasidegenerate. For (iT x ,TT y ) = (i,^) they have 
the same quantum numbers (r,m,n); for (ir x ,ir y ) = (— , — ) and (n x ,ir y ) = (+,+) they differ by 1 
in the quantum number n. 

Upon the transition from OA to OB we have — Q u 3> h giving a u ~ 2 — tt x — ix y . For the 
quantization of the action component J u we choose n = 10 and (ir x ,ir y ) = (— , — ). The transition of 
the WKB lattice takes place in the components J\ and J^, see Fig. |l6|d. In contrast to the transition 
from P to OB the change of the quasidegeneracy is not connected to the action components 



presented in the figure. For OA and OB the same pairs of states are quasidegenerate, see Tab. [V 



Only the shift of the WKB lattice relative to the simple lattice (nh) changes, see Fig. 14. For 
{^x-, %) = (+, — ), (— j +) and (+, +) the pictures are similar and are not shown here. 



Note that we represent the parity boxes in Fig. 16 with side length h/2 also in the direction 
of J^. The reason is that within each region presented in the plots we always have 0a > 0^ or 
©A < 0^> respectively. These relations restrict the range of a^, see Eq. (101), and halve the parity 



boxes P-k. One of these relations is always trivially fulfilled in region O and OB because of the 
different signs of the tunnel integrals there. For OA and P none of these relations holds within the 
whole region, i.e. it is not possible to define some kind of reduced parity boxes that are halve the 
P-n- valid for regions OA and P although this is possible for regions O and OB. 

The situation is much more complicated in the neighbourhood of the intersection line of the 
separatrix surfaces in Fig. ||. Here the transition of the WKB lattice cannot be shown in two 



dimensional sections. In Fig. |17| we present the surfaces X-k v -k z , Hir y -n z and ^7r u7 r z with the action 



component J\, J M or J u , respectively, being semiclassically quantized for the parity combination 
7r = (— , TT y , ir z ) and the corresponding quantum number from (r,m,n). The surfaces carry the 
intersection lines as explained above. We restrict the representation to tt x = — to keep the picture 
clear. For ir x = + the picture is similar. The index * indicates two surfaces with different parities 
located almost on top of each other. With tt x being fixed the effective Maslov phases a\ and a u 
depend only on one further parity ir y or ir z , respectively, depends on ir y and ir z . We therefore 
show two surfaces for J\ and J v and four surfaces for «/„. The parity boxes are omitted in this 



figure. In Fig. 17 thus not every intersection point corresponds to a semiclassical eigenvalue. The 
semiclassical eigenvalues may be identified with those intersection points lying closest to the exact 
quantum eigenvalues represented by spheres in the figure. 



Triaxial Ellipsoidal Quantum Billiards 



30 



a) 



10 
J n 9 
8 



O 



++- 



-+- 



/: 

J \ 



/ 



/ 



-+- 



++- 



-+- 



OA\ 



8 9 10 11 



12 

J, 



13 14 15 16 



b) 



24 

J a 23 
22 



O 



-++. 

Q 



9* 



-++ 



;— + 



36 37 38 39 40 41 42 43 44 

J„ 



c) 



50 



48 



p 


+— 














/ 

/ 


+ 








+— 






++- 

, 4 




" + - 

6 




++-; 

o- 






/ 


r— 
++ 




'-+- 




o 

++- 




-+- 






+ - + 




— + 




o; 

+ + 




- - 




+-+ 

b 




— + 




+-+ 




— + 


+++ 

A 




-++ 




+++ 




++/ 
/ 


t 


+++ 




t>" 
:-++ 




+++ 




i 

-++ 


OB 



24 25 26 27 28 29 30 31 32 



d) 



14 



h 13 



OA 



— + 



^ - + 



12 



-+ 



OB 



8 



9 



10 11 12 13 14 15 16 



Figure 16: Transition of the WKB lattice upon crossing the classical separatrix surfaces (bold dashed lines) in 
action space. In each figure the parity boxes for four neighbouring quantum cells A J\ = A J M = A J„ = 2 are shown. 
The exact quantum states are shown as circles. The semiclassical eigenvalues appear as the intersection of the thin 
short dashed lines within the parity boxes. J is measured in units of h. 
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6 Degenerate Ellipsoids 

So far we treated the billiard in the general triaxial ellipsoid in its classical, quantum mechanical, 
and semiclassical aspects. The triaxial ellipsoid degenerates into simpler systems when any two 
or even all of the semiaxes coincide. In the latter case we obtain the sphere, in the former case 
prolate or oblate ellipsoids which are rotationally symmetric about the longer or shorter semiaxis, 
respectively. In this section we want to take a short look at these degenerate cases where the focus 
is on the similarities in the classical, quantum mechanical, and semiclassical treatment. 

The main theme is that the coalescence of semiaxes of the ellipsoid induces the collision of 
roots or poles. In the classical treatment the disappearence of certain types of motions is expressed 
by the fact that some roots of a hyperelliptic curve collide and the genus of the curve drops. In 
the quantum mechanical treatment it is the singularities of the Helmholtz equation that coalesce, 
which is usually called confluence. Let us look at these transitions in more detail. 

First of all we discuss the ellipsoid itself, without any dynamics. The general ellipsoid with 
parameteres < b < a < 1 has the semiaxes 1 > Vl - b' 2 > Vl -a 2 . There are two choices of 
degenerate cases, either b = a or b = 0. For b = a the longest semiaxes is 1 and the two shorter 
semiaxes coincide, which gives the prolate ellipsoid. In the case 6 = the shorter semiaxis 
is singled out and we obtain an oblate ellipsoid. If a = b = we obtain the sphere. 

The algebraic treatment of the degenerate cases for the classical dynamics is quite simple. In 
the general case the hyperelliptic curves 1Z W have four fixed roots at0<6<a<l and two movable 
roots (i.e. depending on the initial conditions) < s± < a and b < S2 < 1 with s\ < S2- Placing 
the movable roots into the three intervals marked by the fixed roots we obtain four possibilities, 
corresponding to the four types of motion. 

Let us first consider the prolate case where b = a and the ?y-range has vanished. In terms of 
the ranges for £ and £ only case P remains, because the other three become special cases of it, see 



Fig. 11. Only the ordering of the roots < s\ < a 2 < s 2 < 1 is left which means that there is only 
one type of motion in the prolate ellipsoid. The genus of the hyperelliptic curve has dropped to 
one and the remainder of the 77-interval is a pole at z = a 2 , i.e. the integral in Eq. ( |50| ) is elliptic 
and of the third kind. Integrating around this pole gives the residue of this pole divided by 2ir, 
which is 

I v = If, = ±VK 2 -L 2 - 2a 2 E = L x . (107) 

This is the angular momentum of the rotational degree of freedom, which is itself an action because 
the corresponding angle is cyclic. Hence even though the //-interval disappears, the 77-action of 
course does not, this is the reason for the "appearance" of the pole. 

In the oblate case b = the ^-interval vanishes so that in terms of the ranges for £ and 77 the 
two cases OA and OB remain with the corresponding orderings of the roots < s\ < s 2 < a 2 and 
< s 2 < a 2 < s?, < 1, respectively. O becomes a special case of OA and P a special case of OB, 



see Fig. 11. Again the hyperelliptic curve attains a double root, so the genus drops and we obtain 
an elliptic differential of the third kind. The residue at the pole at z = b = gives 2irl u , where I u 
is the angular momentum about the symmetry axis; we find 

I C = I U = L = L z . (108) 

If we finally collapse to the sphere with b = a = 0, there is only the ^-interval left. The motion 
in this interval describes the radial dynamics. The genus of the curve has dropped to 0, because 
now s\ is fixed at zero. The two angular degrees of freedom are hidden in a pole at 0. Considering 



Eq. (13) we see that I = for a = b = so that in this case we find 

I = I v + I lt = K = \L\. (109) 

The fact that there is only one pole while we would like to obtain two actions can be taken as an 
indication that the system is now degenerate, i.e. it has more constants of motion than degrees of 
freedom. 
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Let us now consider the separated Helmholtz equation. The partial fraction decomposition of 
g/ f in Section ^] shows that we have regular singular points at ±a and ±6 in the case of a general 
ellipsoid. The solutions with exponents 1/2 and gave the different parities. In the prolate case we 
have a confluence of a with b and of —a with —b. The equations for £ and £ reduce to the prolate 
spheroidal wave equations. Scaling the variables according to £ = a£ and £ = a£ gives them the 
familiar appearance (see p9]) 



2 

m 



with parameters 



(f _ iy + 2£>' - (A - c 2 f 2 + ~F— = 0, (110) 

r - 1 

(l-CV'-2C>' + (A-c 2 C 2 -^ = (111) 



\={4kE-2a l E)/h% c=2a l E/ff, m 2 = L 2 /h 2 . (112) 



The variable r/ can be turned into an angle after some scaling to compensate for the vanishing of 
the 77-range. The corresponding differential equation yields the familiar result m £ 7L which means 
that each energy eigenvalue is twofold degenerate. Equations ( |110|) and (111) are identical. The 



way of writing them just indicates that they are considered on the different ranges £ G [—1,1] and 
£ > 1. The indicial equations for the regular singular points £ = ±1 and £ = ±1 are of course the 
same and give the exponents ±a = m/2. Note that half the residue of the classical action integral 
over the coalesced singularity gives the exponent of the indicial equation in the quantum case. 

For the oblate case b = we obtain a confluence of b with —b which gives a regular singular 
point at for the equations for £ and n. The two regular singular points at dba usually are removed 
by transforming £ and n separately according to £ 2 = a 2 + a 2 £ 2 and rf = a? — a 2 fj 2 . This gives the 
familiar pair of oblate spheroidal wave equations (see again ) 



2 

m 



(| 2 + 1)^ + 2^'- (A -c 2 ! 2 -^— )V = 0, (113) 

z + 1 

2 

(1 - f?V - 2f)i>' + (A + cV - = (114) 

1 — i] z 

with A and c 2 again defined as in Eq. ( |ll2j ) but now m 2 = L 2 /h 2 . Similarly to the prolate case the 
variable £ can be turned into an angle and the corresponding equation again gives m € Z, i.e. the 
twofold degeneracy of the energy eigenvalues. The indicial equations for the regular singular points 
£ = ±z and rjf = ±1 which correspond to the original regular singular points £ = and n = are 
again the same and again give the exponents a = ±m/2. For the sphere let a = b = which gives 
/ = 0. The scaling £ = rh/^/2E then turns the equation for £ into 

rV' + 2rtp' + (r 2 - n(n + 1))V> = (115) 

with n(n + 1) = \L\ 2 /h?. The variables r\ and £ can be transformed into the azimutal and polar 
angles of the sphere and the corresponding equations give the familiar result n E No and the 
(2n + l)-fold degeneracy of the energy eigenvalues. Eq. (|115|) obviously is the defining equation for 
spherical Bessel functions. It has a regular singular point at with exponents n and — (n + 1). The 
corresponding solutions are the so called spherical Bessel functions of first and second kind. The 
asymptotics at picks out the functions of first kind as the physical solutions for the sphere. 

Note that all the equations cited here for the degenerate cases have an irregular singular point 
at 00. In the non-degenerate case 00 is a regular singular point of Eq. (|56|). 

Concerning the separation of the separation constants the non-degenerate Helmholtz equation 
presents the worst case because it is not separable in the separation constants. In the prolate and 
oblate cases only the constant m can be separated off, the equations are called partially separable 
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for separation constants. The billiard in the sphere belongs to the simplest class of equations which 
are completely separable for separation constants. 

For the semiclassical treatment it is again illuminating to have a look at what happens to the 
hyperelliptic curve TZ W in the degenerate cases. In the non-degenerate case 1Z W has genus 2 and 
therefore has two complex periods giving the penetration integrals 0£ and 0£ in Equations ( j53| ) 
and (|54|). In the prolate and oblate limiting cases one of the handles in Fig. |9] vanishes. The genus 
of the curve drops to 1, i.e. the curve becomes elliptic. Generally an elliptic curve has one complex 
period that gives one penetration integral. In our case it is more useful to think of the elliptic 
curves for the prolate and oblate cases as singular limits of a hyperelliptic curve for the following 
reasons. During the prolate limiting process the middle handle in Fig. |9| shrinks. The penetration 
integrals are not only defined for any b < a but even for a = b where they become infinite. This 
means that although the curve corresponding to the prolate billiard motion is elliptic there is no 
tunneling in our semiclassical treatment. During the oblate limiting process the upper handle in 
Fig. |9| shrinks. Again the penetration integrals are even defined for the limiting case b = where the 
curve becomes elliptic. Here 0^ diverges and 0£ stays finite. The prolate and the oblate limiting 
cases have in common that the penetration integrals connected to the vanishing handle diverge. But 
the oblate billiard still has one finite penetration integral that gives the semiclassical description 
of the tunnelling between tori corresponding to the two types of motions present here. This is why 
the oblate limiting behaviour may be considered as the more typical case. The prolate ellipsoidal 
billiard is peculiar in this sense. This peculiarity is also reflected by the fact that the prolate 



ellipsoidal billiard exhibits quantum monodromy, see [40|. From the point of view of periods of the 
Riemann surface it is important to mention that in the prolate case the sum of the penetration 
integrals is finite. 

The degeneracy of the energy eigenlevels can semiclassically be understood by inspection of the 
energy surfaces. In the prolate and oblate cases the energy surface is symmetric with respect to the 
sign change of the action corresponding to the conserved angular momentum |7]]. Thus the EBK 
quantization condition is each time fulfilled simultaneously at two points of the energy surface. 
The additional degeneracy of the billiard in the sphere classically manifests itself in the resonance 
of the azimutal and polar angular motion. The ratio of the corresponding frequencies has modulus 
1, see Q. The energy surface is thus foliated by straight lines what makes the EBK quantization 
condition being fulfilled not only at one point of the energy surface but on a whole line of points. 
This gives the (2n + l)-fold degeneracy of the energy eigenlevels. 

7 Conclusions and Outlook 

The last section demonstrated the unity of classical, semiclassical and quantum mechanical treat- 
ment in the complex plane, which is one main aspect of this exposition. Another one is to emphasize 
the simple and nice picture of the quantum mechanics of an integrable system as a discretization 
of classical action space. Away from the separatrix surfaces the discretization gives regular lattices 
due to the applicability of the simple EBK quantization of Liouville- Arnold tori in Eq. (0). It is 
much harder to give a semiclassical description of the quantum states whose eigenvalues in action 
space lie close to the separatrix surfaces because of the presence of quantum mechanical tunneling 
between tori with different Maslov indices. The tunneling was incorporated by a uniform WKB 
quantization scheme. This approach is necessary because the application of the simple quantization 
rule ([§) close to the separatrix surfaces in action space can give erroneous additional eigenstates 



[32]. The ingredients for the WKB quantization scheme, i.e. the three classical actions and the two 
penetration integrals, have a consistent interpretation as the real and purely imaginary periods of a 
single Abelian differential of second kind on a hyperelliptic curve of genus 2. In this sense quantum 
mechanics appears as a "complexification" of classical mechanics. 

For the billiard in the ellipsoid we were able to represent all quantum states as a regular WKB 
lattice in the space of the slightly modified actions J which are twice the actions / of the symmetry 
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reduced system. This is impossible in the space of the original actions I. The two classical senses 
of rotations of motion types P, OA, and OB give two different tori in phase space whose actions 
differ in sign. Quantum mechanically we cannot distinguish between these tori and it is therefore 
impossible to assign different quantum state to them. In [1C] all systems like the ellipsoidal billiard 



which (after some symmetry reduction) allow for a representation of the eigenvalues in action space 
are refered to as "one-component systems". In order to represent quantum states in action space 
it is generally necessary to modify the original actions. For systems which are no one-component 
systems the WKB lattice will be much less regular then. 

For future work on ellipsoidal quantum billiards we want to mention two directions. On the 
one hand the ellipsoidal quantum billiard can be taken as the starting point for computations of 
non-integrable quantum billiards resulting from slight distortion of the ellipsoidal boundary. On 
the other hand the semiclassical analysis in terms of periodic orbits is still to be worked out. For a 
chaotic system the Gutzwiller trace formula gives a semiclassical expression for the quantum density 



of states as a summation over isolated periodic orbits [41, 42, 43 1. Analogously the quantum density 
of states of an integrable system can semiclassically be written as a summation over resonant tori, 
i.e. over families of periodic orbits. The Berry- Tabor trace formula gives the quantum density of 



states of an integrable system with / degrees of freedom as a summation over resonant /-tori [44|. 



The main contribution to the density of states stem from the shortest periodic orbits. Generally the 
shortest periodic orbits do not foliate /-tori but lower dimensional tori. Since these contributions 
are not included in the generic case discussed in p4| ] they demand special considerations In 
addition to that the presence of separatrices demands a modification of the Berry- Tabor trace 
formula [32]. Both the non- generic contributions of resonant 2-tori and of isolated periodic orbits, 



and the presence of the crossing separatrices complicate the periodic orbit summation for ellipsoidal 
quantum billiards. For three-dimensional billiards as models for nuclei the shortest periodic orbits 
are important for the explanation of shell structures. For rotationally symmetric billiards this is 
considered in [46, 47 1. For non-symmetric ellipsoids this is still to be done. 
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